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ABSTRACT 

The  usefulness  of  a  barotropic  prediction  model  for  providing  real- 
istic displacements  at  levels  other  than  500  mb  is  suggested  to  circumvent 
difficulties  with  current  baroclinic  models  attempting  to  predict  develop- 
ment as  well. 

A  prediction  equation  is  obtained  by  making  use  of  concepts  from 
Charney's  equivalent  barotropic  model,  'Arnason's  stratified  model,  and 
Cressman's  semi-empirical  barotropic  model.  Further,  a  stream-function 
prediction  scheme,  making  use  of  geostrophic  and  non-divergent  wind  rela- 
tionships, is  used  to  make  a  series  of  24-  and  48-hour  forecasts  of  300-mb 
heights.   The  vorticity  equation  used  is 

^  K>tk)¥t  +  J(V,  Ktt+f)  -Kvffe  J(ij>,V)«  o 

where  K]_,  K2,  and  Ko  are  empirical  constants,,  Results  are  presented  for 
the  effects  of  adjusting  K2  and  Ko. 
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Haltiner  of  the  U„  S.  Naval  Postgraduate  School  for  his  guidance  and  en- 
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-  The  relative  vorticity.   £*  --  jj?  •  ^  x  V. 

ij   -  The  absolute  vorticity.  Tj   --    £  +•  -f . 

■f   -  The  Goriolis  parameter,  Pllsw^,  where  tp   is  the  geographic  lati- 
tude. 

Tm  -  The  mean  value  of  the  Coriolis  parameter,  Q  SL  sw  IS. 

u)     -  The  vertical  wind  component  in  the  x,  y,  p,  t  coordinate  system. 

5t  ' 

q  -  The  upward  component  of  the  apparent  gravitational  acceleration. 

(|5  -  The  geopotential.   <£  =  g  2  . 

'i))  -  The  stream  function  for  the  non-divergent  component  of  velocity. 

■\7  -  The  horizontal  Laplacian  operator  on  a  constant  pressure  surface. 

SJT3-  The  finite-difference  equivalent  of  S7* 

T  (#*>)-  The  horizontal  Jacobian  operator.  J  (A,B)  =  -r~  -r-  —  - —  -rr  . 

c)X  oM    du  OX 

MA,B)-  The  finite-difference  equivalent  of  J  (A,B). 
^J  -  The  finite-difference  equivalent  of  the  derivative,  — —•  , 


1.  Introduction. 

Up  to  the  present  time,  most  attempts  to  predict  the  pressure  height 
at  levels  other  than  500  mb  have  been  made  with  baroclinic  models „  The 
results  have  been  disappointing  because  of  inadequately  controlled  devel- 
opment mechanisms o  Observing  the  persistent  skill  of  barotropic  forecasts 
at  500  mb,  'Arnason  [1]  suggested  the  possiblity  of  devising  a  barotropic 
prediction  model  to  be  applied  to  other  levels,  "capable  of  realistic  dis- 
placements of  pressure  systems,,  „» but  containing  no  mechanism  for  develop- 
ment . " 

This  investigation  introduces  an  empirically  adjusted  divergent  baro- 
tropic model  for  application  to  the  300-mb  surface.  Concepts  from  Charney's 
[2]  equivalent  barotropic  model  and  ' Arnason' s  stratified  model  [1]  are 
incorporated  in  addition  to  empirical  results  by  Cressman  [33  on  control 
of  very  long  atmospheric  waves. 


2.  Background. 

Looking  first  to  the  equivalent  barotropic  model,  it  is  recalled  that 
partial  baroclinicity  is  permitted  in  that  wind  speed,  but  not  direction, 
is  allowed  to  vary  with  height  at  any  point „  If  the  integrated  mean  of 
this  wind  with  respect  to  pressure  is  determined,  the  wind,  V ,  and  vor- 
ticity,  g   ,  at  any  level  can  be  expressed  as  a  scalar  factor,  A,  of  this 
mean.  For  further  simplicity,  A  is  considered  to  be  a  function  of  pressure 
alone .   Thus 

V  =  Acp)V  ,   X   =  Acp>?  ,  (2.1) 

where  the  cap  indicates  the  vertically  integrated  mean.  Climatologically, 
^V  is  observed  to  occur  at  about  600  mb0 

Substituting  the  expressions  (20l)  in  the  approximate  vorticity 
equation 

ft  *v.v(jr*o«k£ 

and  integrating  with  respect  to  pressure  gives 
/s 
||  4-  V-  V(Adr  +  0  =0.  (2.2) 

Here  f  is  the  coriolis  parameter,  and  A  is  found  empirically  to  have  a 
value  of  about  1,25°  Multiplying  (2.2)  by  A  and  substituting  A  *f=  £  and 
A  V=V  forms  the  basis  for  barotropic  forecasts  at  the  equivalent  baro- 
tropic level »  (p";;")i)  occurring  in  nature  at  about  500  mb0 

If,  however,  the  expressions  (2.1)  are  substituted  for  V  and  £  into 
(2.2),  the  following  is  obtained: 

|f  +  V-^(Ka?+f)  =0   ,  (2.3) 


where  K<>  =  -—,  For  the  wind  field  (2,1),  equation  (2,3)  applies  to  any 

A 

level  with  the  Ko  appropriate  to  that  level. 

In  applying  the  theory  of  the  stratified  model  to  forecasting, 
1 Arnason  arrives  at  an  expression  for  divergence 

*■* --MS +  *-**)  >  (z'h) 

where  <B  is  geopotential,  V  is  the  basic  flow,  and  H  is  the  thickness  of 
the  stratified  layer.  The  parameter,  g1 ,  is  related  to  gravity,  although 
much  smaller,  and  depends  upon  static  stability  and  the  vertical  extent 
of  the  stratified  layer.  This  expression  for  divergence  is  introduced  into 
the  vorticity  equation 

to  give 

ff  +  V-Vr,,  _2L(H  +  v.V$)=0,  (2"6) 

where   Yj  =^++      is  absolute  vorticity;  "Y)a  is  defined  for  our  purposes  as 
Y]=   KA?  I   f       of  (2,3). 

The  prediction  scheme  consists  of  introducing  a  stream  function,  "W  , 
to  represent  the  wind  and  vorticity 

V=l?xVV    ,         X   =  Va  V  (2-7) 

where  lr  is  the  unit  vector  in  the  vertical0  The  stream  function  is  ini- 
tially obtained  from  the  balance  equation 

£vay  +5(v«vaj  -  v,y)  +  vy-vf  -  va$=o  .  (2.8) 

Next,  in  order  to  simplify  the  procedure,  the  quasi-geostrophic  relation- 


ship 


$  =  fV  (2.9) 

is  introduced  into  the  divergence  term.  The  result  of  substituting  (2„7) 
and  (2,9)  into  (2,6)  is 


v*y  t  J<v,>)a)-M$  -  J(V,V)]  =o 


or 


^-f^)S+J(V,^)  "^  J(^)-0,  (2-10) 

where 'y  is  the  stream  function  corresponding  to  the  basic  flow.  Equation 
(2.10)  is  a  prognostic  equation  in  a  single  variable, ^  . 


3 .  Procedure , 

Application  of  'Arnason's  stratified  model  involves  an  interpretation 
of  the  parameters,  V  ,  g9 ,  and  H.  Also,  the  introduction  of  a  factor,  K]_, 
into  the  tendency  portion  of  the  divergence  term,  and  a  factor,  Ko,  into 
the  advection  portion  of  the  divergence  term  permits  adjustment  of  the  mag- 
nitude of  these  terms  for  various  levels. 

The  basic  flow  for  this  prediction  model  is  taken  as  a  spaced-averaged 
flow  at  300  mb,  employing  a  15-pass,  Appoint  smoother  over  the  finite-dif- 
ference grid  used  for  computations  (grid  distance  equal  to  381  km  at  60N 
latitude): 

Yo  -  7J-(V.  +  U  +  Va  +Vh)  (3a) 

The  subscript  zero  denotes  the  central  point  of  a  diamond  array  and  the 
subscripts  one  through  four  denote  the  surrounding  four  grid  points. 

In  arriving  at  an  interpretation  of  g' ,  reference  is  made  to  Cressman' s 
[33  treatment  of  a  model  including  a  dependence  upon  static  stability.  He 
introduced  a  stability  parameter,  a.  ,  which  as  applied  here  amounts  to  a 
substitution  of 

t  -   F  (3-2) 

In  a  series  of  500-mb  forecasts  Cressman  used  various  values  of  zx.  in  his 
divergent  barotropic  model*  It  was  concluded  that  P-   —   4  affords  the  de- 
sired control  on  the  ultra-long  waves,  and  gives  better  forecasts  than  a 
former  Joint  Numerical  Weather  Prediction  Unit  forecast  model.  Use  has 
been  made  of  this  value  by  the  substitution  of  (3»2)  into  (2o10)« 

In  addition,  the  mean  value  4? fSSP has  been  used  for  -S-P-  (similar  to 
Cressman),  taking  z5oo  as  the  U„  Sc  standard  atmosphere  500-mb  height, 


and  fm  =  f  at  45N  latitude. 


With  these  last  adjustments,  (2,10)  becomes 

(V*_K|^)^  tj(v,K,r+o-K,^»jfi'y)-o.     (3.3) 

In  finite  difference  form,  (3°3)  becomes 

&-K^€)ft   +  4  9CV/A0  -  Ky^g J(?,V)=°  ,   (3-4) 

where  SFa,  A  ,  and  Q  are  the  finite -difference  equivalents  of  the  Lapla- 
cian,  differential,  and  Jacobian  operators,  respectively;  and  $  is  the 
sea-level  distance  between  grid  points .  Forecasts  are  made  in  one-hour 
time  steps,  using  a  forward  difference  for  the  first  step  and  centered 
differences  for  successive  steps. 

Equation  (3.4)  was  programmed  for  the  Control  Data  Corporation  Model 
1604  computer  using  U.  S.  Fleet  Numerical  Weather  Facility's  1977-point 
octagonal  grid,  which  covers  the  Northern  Hemisphere  from  pole  to  about 
ION. 

Step  one  in  the  forecasting  procedure  consists  of  solving  the  balance 
equation  for  an  initial  stream  function,  using  analyzed  D-values  (deviation 
in  height  from  standard) „  A  technique  formulated  and  programmed  for  the 
CDC  1604  by  "Arnason  and  Reese  at  the  Fleet  Numerical  Weather  Facility  was 
used.  It  was  found  that  for  most  of  the  300-mb  data  available  for  this 
investigation,  this  method  would  not  converge  to  the  cut-off  criterion 
(all  residuals  less  than  0.5  feet  in  the  relaxation  procedure  as  normally 
used  at  lower  levels).   In  order  to  proceed  with  the  work,  an  arbitrary 
cut-off  was  introduced  if  non-convergence  became  apparent.  This  intro- 
duced errors  at  a  few  points  of  the  grid  due  to  residuals  which  were  now 
of  the  order  of  three  to  five  feet.  These  were  considered  to  be  insignif- 


icant  with  respect  to  the  prediction  scheme. 

Step  two  was  solution  of  equation  (3«4),  and  this  was  repeated  on 
each  day's  data  using  various  values  of  the  factors  K2  and  Ko„   The  factor 
Kt  was  taken  equal  to  one  and  not  adjusted  in  this  investigation .  Values 
of  K2  and  Ko  were  chosen  as  work  progressed  so  as  to  reduce  forecast  er- 
rors to  a  minimum. 

Step  three  consisted  of  verifying  the  forecasts  against  the  prepared 
analyses  by  a  simple  point -by-point  difference,  and  calculation  of  the 
pillow  and  root-mean-square  error  (RMSE): 

P.llow  =i^(A-8)n   ,  RMSE  =[i£(A-B-  Pillow)*"       (3.5) 


H.-I 


H=l 


x  =  1977  points,  A  =  forecast  map,  B  =  analysis  map, 


4.  Results. 

Table  4.1  lists  the  results  of  adjusting  the  factor  K3  with  K2  held 
equal  to  one.  Numbers  to  the  left  of  each  column  are  pillow  while  the 
numbers  to  the  right  are  RKSE.  Upper  numbers  of  each  paired  set  are  for 
24-hour  forecasts  and  lower  numbers  are  for  48-hour  forecasts „ 

Runs  were  not  made  for  all  included  values  of  K3  but  minimum  errors 
were  sought.  The  numbers  in  parentheses  are  values  taken  from  faired 
curves  through  the  values  that  were  run.   Fig.  4.1  shows  the  curves  for 
24-hour  forecasts j  those  for  48-hour  forecasts  would  be  similar .  The 
trends  of  errors  for  the  successive  values  of  K3  were  observed  to  be  quite 
regular,  thus  justifying  an  approximation  of  intermediate  values  where 
necessary  for  the  averages  over  all  days.   In  cases  where  actual  reversal 
of  error  trends  was  not  observed,  the  fact  that  differences  between  suc- 
cessive values  of  K3  were  reduced  to  zero  or  one  foot  was  taken  as  an  in- 
dication that  reversal  was  imminent,  and  ad justments  of  K3  were  terminated. 

In  one  case  (17  January)  a  forecast  was  made  using  a  positive  value 
of  K3,  showing  an  error  much  larger  than  the  optimum  for  this  day.  Other 
work  on  this  term  not  presented  here  further  confirms  the  large  errors  for 
values  of  K3,  including  K3  equal  zero  (exclusion  of  the  term),  beyond  the 
range  indicated  herein. 

The  column  entitled  "Balancing  errors"  indicates  the  differences  be- 
tween (a)  the  height  Tields  obtained  after  balancing  to  the  stream  func- 
tion and  inverting  again  to  height,  and  (b)  the  initial  height  field. 

Table  4.2  lists  the  results  of  adjusting  the  value  of  I^?.  Since  K3 
was  first  adjusted,  optimum  values  of  this  factor  were  used  in  adjusting 


Kp  • 


It  will  be  noted  from  table  4.1  that  a  considerable  portion  of  the 
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forecast  errors  is  apparently  due  to  the  stream-function  conversion.  On 
13  January,  the  conversion  error  amounts  to  Ul%   of  the  total  forecast  er- 
ror, while  on  other  days  the  ratios  of  errors  range  from  17%   to  25^. 

The  trends  of  errors  for  varying  K^  all  indicate  a  minimum  error  in 
the  range  -0.5  to  -1.75  and  also  indicate  the  advantage  of  including 
'Arnason's  divergence  term.   From  extrapolation  of  the  curves  in  fig.  4.1, 
errors  for  Ko  equal  zero  are  estimated  to  average  about  20  feet  above  the 
least  average  for  24-hour  forecasts. 

The  use  of  the  factor  K2  does  not  seem  to  justify  any  deviation  from 
unity,  with  improvements  ranging  from  zero  to  only  three  feet  for  24-hour 
forecasts c 


CTTD 

ITN 

Cxi    C^ 

O  rH 

H  t>- 

nO  CO 

NO   On 

m  in 

t>- 

CO  O 

x£>    UA 

IT\  r>- 

r>-  CV 

IT\  H 

nO  o 

• 

CV  -4 

CNJ    C>-\ 

cv  en 

CV  -4 

CV  -4 

cv  en 

H 
1 

O  en 
CV  CV 
+    + 

O  H 
CV 

+    1 

O  -4 

•          • 

o  m 

H  en 

CO  CO 

V\-4- 

rH  CV 

On  -4 

m 

CO  o 

vO    m 

en  in 

nO  o 

ir\  o 

UTN.C0 

• 

cv  -4 

cv  en 

cv  en 

CV  -4 

cv  -4 

cv  en 

H 

N^->"*— * 

1 

0--  en 

O  H 

O  H 

O  -4 

CV  CV 

cv 

-4  -4- 

H  CV 

+  + 

+   1 

+  + 

+    + 

O  nO 

m 

o  cv 

en  o 

'S  ncT 

CO  H 

O  -4 

nO  CO 

CM 

CO   r-H 

sO  o 

en  -4 

V\  On 

-4  On 

U"\  f- 

* 

cv  -4 

cv  en 

cv  en 

cv  en 

cv  en 

cv  en 

H 

n. y*-^s 

1 

On  m 

O  H 

en  o 

O  ~4 

cv  cv 

cv 

iHnO 

rH   CV 

+  + 

+   I 

+  + 

+    + 

CO  nO 

P  \  cv 

CO  H 

U>  CO 

ut\nO 

CO  <H 

ir\  H 

O 

CO  cv 

nO  c- 

cv  en 

into 

-4  On 

in  co 

H 

cv  -4 

cv  en 

CV  en 

cv  en 

cv  en 

cv  en 

O  c*A 

o  cv 

o  cv 

^Q 

O  -4 

cv  cv 

CV 

UN.    -4 

HnO 

iH  CV 

+  + 

+    1 

+  + 

+  + 

+   + 

en 

W 

COnO 

U"\ 

C-~vO 

irwo 

CV  u"\ 

l/N,  IT\ 

O  H 

C^nO 

r- 

CO  en 

f-  CO 

cv  en 

IT\  CO 

U"\  O- 

if\  CO 

« 

cv  -4 

cv  en 

CV  en 

cv  en 

cv  en 

cv  en 

o 

O  en 
u~\  -4- 

en  On 
H  wr\ 

O  -4 
H  CV 

+  + 

+  + 

+  + 

CV  nO 

•        • 

CM  cv 

-J"  U"N 

H  O- 

CO  no 

nO  en 

CV  -4 

ITS 

O  ir\ 

CO  o 

CV  en 

uaCO 

in  o 

nO  On 

• 

cv  -4 

cv  -4 

CV  en 

cv  en 

CV  en 

CV  en 

o 
1 

CO  in 
cv  cv 

+  + 

H   CV 
CV 

+    1 

O  -4 

+  + 

ir\ 

nO  -4 

CV 

H  On 

en  -4 

rH 

+ 

-4  CO 
H  ir\ 

+  + 

bO 

.5  g 

°  & 

o 

cv 

H 

en 

o 

en 

c  p 

u-\ 

U"\ 

O 

nO 

m 

nO 

rt   H 

"-*    r^! 

o 

■4 

ON 

nO 

cv 

CI)    W 

H 

en 

cv 

<H 

CO 

+ 

+ 

+ 

+ 

+ 

CO 

CD 

o 

o 

c 

c 

c 

t£ 

CD 

CD  rH 

CD  H 

trj  cv 

crj  cv 

rd  CV 

0} 

-P 

a  no 

Q  nO 

-D  n0 

•"3  nO 

•s  nO 

u 

0} 

o 

ON 

On 

On 

ON 

CD 

Q 

-4-  H 

ir\  H 

en  H 

r-H 

H  rH 

> 

H 

rH 

•H 

rH 

CV 

< 

CD 
rH 

CCJ 


cn 
c 

•H 
-P 

DO 

-a 
o 

CO 

-p 
H 

CO 

CD 


10 


\ 


\ 
\ 
\  \ 


\  \ 


320 


3(0 


300 


,l7Jan 


\ 


\      W 


1  I 1 I I 


220 


+1.0  +  .5  0  -.5 

K3 

Fig,   4.1-Trencls  of  RMSS  values,   24-hour  forecasts 

11 


-1.0  -Lb' 


i^ 

t^- 

O-  CO 

o    c*--4 

CO  CO 

vO  o 

r--- 

cv  en 

II     cv  C^ 

O 

o  -4 
en  en 

+   + 

cv  en  rH 

—     +    + 

O  en 

nO  CO 

CO   "A 

rH   CO 

cv  <n 

CO  CO 

vO   lA 

CV  C^ 

vO  CC 

LT\  O 

co 

cv  en 

cv  en 

cv  en 

cv  en 

CV  -4 

o 

O  H 

CV   rH 

H  sO 

en  rH 

rH  en 

en  en 

CV 

UA  -4 

rH  vO 

rH  CV 

+  + 

+  + 

+   + 

+  + 

+  + 

cv 

ixl 

t—  o 

o  o 

rH  H 

tf>  -4 

vO  o 

c-  CO 

u~\  -4 

cv  en 

U-NCO 

-4"  o 

o 

cv  en 

CV   en 

cv  en 

cv  en 

cv  en 

r> 

o  r- 

r-1   O 

H  U> 

en  rH 

o  <r 

CV  CV 

CV 

UA  -4 

HMD 

rH  CV 

+  + 

+ 

+    + 

+  + 

+  + 

O  "A 

O  rH 

rH  C*- 

u-\0 

CO    rH 

CO  o 

v£)   UA 

CV  en 

tr\  CO 

-4  O 

CV  -4 

CV   <n 

CV  en 

cv  en 

cv  en 

rH 

O  en 

Or* 

o  -4 

en  o 

O  -4 

CV  CV 

CV 

UA^j- 

HsO 

<h  cv 

+  + 

+   1 

+  + 

+  + 

+  + 

ir\ 

CO 

tf\ 

t> 

ir\. 

o 

o 

— 

♦ 

• 

• 

• 

0 

H 

1 

1 

o 
1 

rH 

1 

rH 

1 

o 

O 

c 

c 

§  cv 

<D 

<D  H 

G)  rH 

05  CV 

efl  CV 

-p 

Q  vO 

Q  vO 

►-3  vO 

-a-O 

•-aso 

o3 

O 

On 

o 

O 

ON 

a 

H*H 

iTSrH 

enrH 

C--rH 

H  rH 

rH 

H 

rH 

cv 

cv 

-4 

0> 
rH 
43 

nj 
Eh 


CV 

tri 

M 

C 
•H 
+3 

CO 

•r"3 

-o 

crj 

<H 

O 

(0 

•P 
H 
P 
(0 

<u 


12 


Figs.  4.2  through  4.4  depict  the  13  January  analyzed  height  field, 
the  14  January  analyzed  height  field,  and  the  13  January,  24-hour  forecast 
height  field  verifying  on  14  January,  respectively.  In  fig.  4.4  the  ini- 
tial (13  January)  and  observed  (14  January)  characteristics  of  certain  of 
the  systems  are  superimposed  upon  the  forecast  (for  14  January)  field  to 
indicate  contrasting  features.  Initial  features  are  dotted,  and  observed 
are  dashed. 

The  following  characteristics  may  be  noted; 

(1)  In  the  western  Pacific,  the  forecast  of  the  short-wave  trough 
was  very  much  in  agreement  with  its  actual  movement  and  intensity.  No 
cyclogenesis  occurred  in  this  system  and  forecast  errors  were  near  zero 
in  the  area  of  the  trough. 

(2)  The  eastern  Pacific  low  occurring  on  13  January  almost  completely 
diss  .pated  in  24  hcurs.  The  forecast  position  of  this  low  is  nearly  that 
of  the  14  January  trough  identifying  the  system.   Height  errors  in  excess 
of  700  feet  are  due  to  cyclolysis,  which  of  course  is  not  provided  for  in 
the  model. 

(3)  In  the  Alaskan  ridge,  cyclogenesis  occurred  from  13  to  14  Janu- 
ary in  the  NW  portion,  accounting  for  errors  up  to  600  feet,   The  ridge, 
however,  was  well  forecast, 

(4)  The  trough  off  the  Atlantic  coast  of  the  U,  S,  was  forecast  to 
move  in  agreement  with  observation,  but  cyclogenesis  again  resulted  in  er- 
rors of  about,  400  feet. 

(5)  The  low  initially  over  the  central  Mediterranean  moved  into  a 
cut-off  position  south  of  a  blocking  high.  The  forecast  movement  was  in 
fair  agreement,  although  a  bit  short,  but  accompanying  cyclogenesis  raised 
errors  to  about  600  feet.  The  forecast  position  of  the  blocking  high  was 
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good,  but  again,  anticyclogenesis  in  the  high  and  the  ridge  to  the  north- 
east introduced  forecast  errors  of  400  feet  and  600  feet,  respectively. 

Elsewhere  on  the  map,  the  significant  errors  seem  generally  to  be  due 
to  development  processes  and  not  to  movements  of  the  systems.  These  ob- 
servations were  corroborated  by  investigation  of  another  sequence  of  maps 
produced  in  the  investigation. 
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5.  Conclusions. 

It  is  noted  from  the  averaged  errors  in  table  4»1  that  there  are  only 
limited  differences  in  results  over  the  range  of  Ko  values  indicated „   It 
was  apparent,  however,  that  sizable  differences  occur  as  Ko  is  extended 
beyond  this  range,  and  of  particular  interest  is  the  advantage  over  exclu- 
sion of  the  term  (Ko  equals  zero) .  These  results  indicate  an  optimum  val- 
ue for  Ko  of  about  -1,0. 

Even  though  varying  K2  did  not  yield  significant  differences  in  this 
investigation,  it  is  possible  that  a  greater  number  of  cases  might  do  so, 
particularly  if  some  other  season  is  considered,,  Only  data  from  December, 
1961  and  January,  1962  were  used. 

Taking  note  of  the  large  errors  encountered  in  converting  to  a  stream 
function,  it  is  possible  that  a  geostrophic  model,  not  requiring  such  con- 
versions, might  reduce  total  errors  somewhat.   It  is  hoped  that  this  in- 
vestigation can  be  accomplished  in  the  near  future. 

The  model  developed  in  this  paper  seems  to  fulfill  its  intended  pur- 
pose, that  is,  forecasting  realistic  displacements  unaffected  by  uncon- 
trolled developments.   *.:ith  an  understanding  of  its  limitations,  a  fore- 
caster should  have  here  a  useful  tool,  when  supplemented  by  techniques  for 
development  prediction,  to  forecast  pressure  systems  at  the  300-mb  level 
more  easily  and  accurately. 
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